Local Behavior of Sparse Analysis Regularization: Applications to Risk Estimation

In this paper, we aim at recovering an unknown signal x0 from noisy measurements y = Φx0 +w, where Φ is an ill-conditioned or singular linear operator and w accounts for some noise. To regularize such an ill-posed inverse problem, we impose an analysis sparsity prior. More precisely, the recovery is cast as a convex optimization program where the objective is the sum of a quadratic data fidelity term and a regularization term formed of the `-norm of the correlations between the sought after signal and atoms in a given (generally overcomplete) dictionary. The `-sparsity analysis prior is weighted by a regularization parameter λ > 0. In this paper, we prove that any minimizers of this problem is a piecewise-affine function of the observations y and the regularization parameter λ. As a byproduct, we exploit these properties to get an objectively guided choice of λ. In particular, we develop an extension of the Generalized Stein Unbiased Risk Estimator (GSURE) and show that it is an unbiased and reliable estimator of an appropriately defined risk. The latter encompasses special cases such as the prediction risk, the projection risk and the estimation risk. We apply these risk estimators to the special case of `-sparsity analysis regularization. We also discuss implementation issues and propose fast algorithms to solve the ` analysis minimization problem and to compute the associated GSURE. We finally illustrate the applicability of our framework to parameter(s) selection on Email addresses: [email protected] (Samuel Vaiter), [email protected] (Charles-Alban Deledalle), [email protected] (Gabriel Peyré), [email protected] (Charles Dossal), [email protected] (Jalal Fadili) Preprint submitted to Applied and Computational Harmonic Analysis October 10, 2012 ha l-0 06 87 75 1, v er si on 2 10 O ct 2 01 2 several imaging problems.